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A DIFFUSION APPROXIMATION MODEL 
FOR A COMMUNICATION SYSTEM 
ALLOWING MESSAGE INTERFERENCE 

Donald P. Gaver 
Naval Postgraduate School 
John P. Lehoczky 
Carnegie-Mel Ion University 

1 . Introduction and Problem Statement 

We study the operating characteristics of an element of a complex commun- 
ication system, the element consisting of c channels which service an 
arrival stream of messages. When a message arrives, it selects a channel at 
random and initiates a transmission (service) time of random duration. If by 
chance the channel selected is already occupied, i.e. is in the process of 
transmission, both messages may be "destroyed," i.e. terminated before com- 
pletion, and the channel reverts to an empty or open condition. The trans- 
mitters of the messages are capable of detecting the event of destruction, 
and following such an event go into a re-try or re-transmi ssion population, 
from which they later make attempts to occupy a channel and eventually com- 
plete message transmission. We develop an approximate probability model to 
describe the performance of such a system. The approximation tends to become 
exact as the number of channels, c , becomes large (c-x») , but numerical 
studies indicate that it may be quite adequate for c near ten. 

The above problem bears close resemblance to a problem of packet switching 
of data on the ARPANET and, very likely, on other satellite communications 
systems, as discussed in Kleinrock (1976), p. 362 ff. We choose to represent 
the various interacting populations of messages, i.e. those in process, and 
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those in the re-try population, by means of stochastic differential equa- 
tions as was done by Gaver and Lehoczky (1976) for a related problem. 
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2. Mathematical Formulations: Model 1 



Assume that messages arrive at a system of c channels according to a 
non-homogeneous Poisson process of rate cA(t) ; that is, the probability 
that a message arrives in (t, t+dt) is cX(t)dt + o(dt) . Each message 
then selects a channel at random; any particular channel , whether occupied or 
not, is selected with probability c"^ . If the channel is unoccupied the 
message begins transmission; the probability that it completes in time 
(t, t+dt) is y(t)dt + o(dt) . Note that if y(t) = y > constant, the 
messages enjoy exponential service times, and if we wish to generate other, 
say Erlang, service times the device of phases, or extra artificial compart- 
ments, may be used, as in Gaver and Lehoczky (1976). 

If a message in progress on a channel is interrupted by the arrival of 
another message, both are assumed instantly destroyed, and the message initi- 
ators are transferred to a re-try population , R . A message that has entered 
R changes its status at time t with probability v(t)dt + o(dt) : with 
probability a(t) the change of status implies an attempt to re-transmit 
on a channel, and with probability a(t) = l-a(t) the change of status 
implies loss -- the message may no longer be worth transmitting. 

Our representation of the above setup is in terms of the following state 
variabl es . 

Q(t) = the number of messages being transmitted, i.e. occupying channels, 

~ at time t ; clearly 0 s _g(t) < c . (2.1) 

^(t) = the number of messages in R at time t ; 0 ^ £(t) < «> . 

L^(t) = the total number of messages that have been lost by time t ; 

~ 0 < L(t) < . 

The state of the system is thus (£(t), R(t), £(t)) , a discrete vector- 
valued Markov process, according to our earlier assumptions. We shall 
characterize this process approximately when c^ , and shall in particular 
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treat £ , ^ , and I. as continuous stochastic processes, in fact as 

diffusion processes, see Arnold (1974). 

Here is a formalization of the transitions described earlier: 



Trans i tion 


Probability 


(in t, t+dt) 

(a. £. L) - (2+1 . R, L) 


a(t) 

cX(t)[l-— ]dt 


+ (a-1. £+2. L) 


2,(t) 

cA(t) dt 


+ (Q-l, R, L) 


y(t)£(t)dt 


* (2+> . £-1 . £) 


Q 

v(t)a(t)[l-J]R(t)dt 


+ (2-1 , £+1 , L) 


a 

v(t)a(t)^(t)dt 


+ (2. £-1 . ui ) 


v(t)a(t)R(t)dt 



Other transitions have negligible probability of occurring in (t, t+dt) . 
Now in principle one can write down Kolmogorov equations for the transition 
probabilities of the (£, l^) process and solve them. However, such a 

solution must inevitably be numerical. Here we shall write down an approxi- 
mate system of Ito stochastic differential equations, see Arnold (1974) > to 
describe process evolution, and from them deduce certain useful information 
valid when c is large. We write, on the basis of (2.2), 

2(t) £(t) £(t) 

dQ(t) = {cX(t)[l-'^] - cX(t)^- p(t)2(t) + v(t)a(t)[l-'^]R(t) 

£(t) 

- v(t)a(t)-^ — £(t)}dt 



wr 



ITtT 



cA(t)[l --^ — 1 dWi(t) - /c^(t)-^ — 



wr 



/y(t)£(t) dW 3 (t) + /v(t)a(t)[l-'^^]R(t) dl^(t) 

/' 

yv(t)a(t)— R(t) dWg(t) , 



(2.3) 
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Q(t) Q(t) 

dR(t) = {2cX(t)^7 v(t)a(t)[l- ]R(t) 

Q(t) 

+ v(t)a(t)^^ R(t) - v(t)a{t)R(t)}dt 

/ ^ ^ 

+ 2/cA{t)^^ dW2(^) - /v(t)a(t)[l- ]R(t) d!;^(t) 



/ Q(t) / 

+ /v(t)a(t)^^ R(t) dWg(t) - /v(t)a(t)R(t) dWg(t) , 



dL(t) = v(t)a(t)R{t)dt + /v(t)a(t)R{t) dW^(t) 



where (W^. (t), o<t, i=l, 2, .. 6) is a 6 dimensional standard Wiener process 
whose components are independent. 

The rationale for writing (2.3) is as follows. Each term in the brack- 

etted part of the expression for, say, dQ , is a component of the drift 

or expected change in Q between t and t+dt ; e.g. cA(t)[l- is 

0 ^ c 

the expected increase in Q caused by a newly arrived message immediately 

Q(t) 

reaching an empty channel (line 1 of (2.2)), while cX{t )'^ — represents 

the decrease in Q caused by a new message that arrives and chooses a busy 

channel, only to be immediately transferred to R , along with the message 

in progress (see line 2 of (2.2)). The coefficients of the Wiener process 

terms are proportional to the standard deviations of the motions in 

(t, t+dt) corresponding to the drift terms. 

Next, introduce an expansion of (Q, R, L) to terms of order / c . 

0*^ 

Define the stochastic noise processes by 
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X(t) 



Y(t) = 



Z(t) = 



Q(t) - cq(t) 



,/c 

R(t) - cr(t) 
j^(t) - c5.(t) 



(2.4) 



the vector (q(t), r(t), l{t)) = lim (£(t)/c, £(t)/c, JL(t)/c) ; the exis- 

C->oo 

tence of these limiting functions is suggested by results of Kurtz (1971). 

An application of Ito's lemma, Arnold (1974), p. 90, provides a sto- 
chastic differential equation (s.d.e.) for ()^, Y^, 1) as follows (hereafter 
we do not explicitly express the t-dependence) : 

X 
Y 
1 



/• \ 

dX 




r 

- (2X + y + 2avr) 


av(l -2q) 






dY 


= 


2X + 2avr 


-av(l-2q) - av 


dZ 




0 

L 


av 



- ^ 



q' - X(l-2q) + yq - avr(l-2q)' 
r' - 2Xq + avr(l-2q) - avr 
I ' - avr 



(2.5) 



/X(l -q) -v^ /avr(l -q) -/avrq 

0 2\/Xq 0 /av(l -q) /vaqr 

0 0 0 0 0 



0 

Vvar 

/var 



~3 

Ha 

ik 

J!6 



e 



Now th6 coefficiGnt of /c must be identically zero, for otherwise our 
system of equations does not converge. This leads to the 
Deterministic Equations : 

q'(t) = X(l-2q) - yq + avr(l-2q) 

r'(t) = 2Xq - avr(l-2q) - avr (2.6) 

' (t) = avr 

The solutions to these equations, which must be obtained by numerical inte- 
gration, provide a deterministic approximation to (Q, R, L) -- in fact, 
(cq, cr, c£) :::: (e[q], e[r], e[l]). 



Stochastic Equations : 

In view of (2.6), (2.5) now appears in the form 



dt + B.dW 



dX 




f \ 

X 


dY 




Y 




L 


0^ 


dZ 




1 









(2.7) 



(a 3 X 2 matrix) being identified as the coefficient of {1, D' ’ 
and (a 3 x 6 matrix) identified as the coefficient of the Wiener process 

term, in (2.5). Next note that the s.d.e. has a special form if the Z 



term is omitted; it is reasonable to focus on X and Y , for clearly 

^-(t)-^®® as t-x» , since represents cumulative losses. Thus we shall be 

interested in the first two equations in (2.6), and in (2.7) which we now 
write as 



dX 




X 




= £t 




dY 




Y 


V J 







dt + D^.dW 



( 2 . 8 ) 
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and C. (2 x 2) is A. without its last row, while D 4 . (2 x 6 ) is B^. 

~t ~t o-t ~t 

without i ts last row. Now the bivariate stochastic process described by 

(2.8) is nonstationary Ornstein Uhlenbeck, since C. has eigenvalues with 

non-negative real parts; see Arnold (1974), Sec. 8.3, for an account of the 

scalar case, or see Schach and McNeil (1973). Much that is useful is known 

about this process. In particular, if (X(0), Y(0)) = 0 , and q(0), r(0) 

are given, then for all t>0 (X(t), Y(t)) has a bivariate normal distri- 

bution with mean 0 and covariance matrix E. which satisfies the differ- 

ential equation 

L = C..E. + E.C; + D.D! (2.9) 

see Arnold (1974), Sec. 8.2. 

Combining these facts leads to the 

Result : (£, R) is approximately bivariate normal (Gaussian), as c-^ : 

(£(t) , R(t)) ~ W(c(q,r) , cE^) (2.10) 

From this expression it is possible to estimate the probability that at 
least a specified number of channels are occupied at time t , and also to 
estimate the probability that there are no more than any specified number of 
customers awaiting re-try. Such quantities are useful measures of system 
performance. Of course both the deterministic equations (2.6) and the 
equation (2.9) for the covariance matrix E. must be solved numerically, 
but this should be a less difficult step than is the simulation of such a 
system. 

Steady-State Results : 

Suppose X(t)->X , y(t)-^y , v(t)-^v , a(t)->a all positive constants, 

as t-w . In this case it may happen that q(t)-»-q and r{t)^r , the latter 
satisfying ( 2 . 6 ) with derivatives set equal to zero: 
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0 = A(1-2q) - pq + avr(1-2q) 
0 = 2Aq - avr(l-2q) - avr 
and these are immediately solved to give 



( 2 . 11 ) 

( 2 . 12 ) 



^ _ (l+2p) - /(l+2p)^- 8ap _ 
^ 4a " 



p-q 
a n 



2p 



(l+2p) + /(l+2p)^- 8pa 



(2.1.3) 



where p = A/y and n = v/p » and provided that the numerical results ob- 
tained are feasible: 0<q<l , 0<r<«> . From (2.10), it is necessary that 

q < ' 2 - in order for a steady state to exist. 

When steady state conditions hold Z.=0 , so the steady state co- 



variance matrix, I , is the unique positive definite solution of 

-DO' = CE + ZC 



(2.14) 



where 



-(2A + p + 2avr) av(l-2q) 

2A + 2avr -av(l-2q) - av 



(2.15) 



and 



DD' = 



A + pq + avr 
-(2Aq + avr) 



writing 



E = 



•(2Aq + avr) 
4Aq + vr 



o 0 

11 12 

o o 

12 2 2 



(2.16) 



(2.17) 



we may express 





r 


a 




1 1 




a 


1 2 




a 

22j 





b(a+b-d) 2b^ b^ 

bd 2ab ab 



2ad a(a+b) - bd 



A + pq + avr 
-(2Aq + avr) 

4Aq + vr 



(2.18) 
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where 



a = 2X + M + 2avr , d = 2X + 2avr 

(2.19) 

b = av(l-2q) , A = 2ab(d-a-b) + 2b^d . 

The values of q and r are available from (2.13). From these results 
actual numerical expressions for the steady state probability distributions 
of Q and R may be computed. Notice that the long-run loss or defection 
rate of the final equation of (2.6) i.e. cl' i°°) ~ cavr , where r comes 
from (2.13). 
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3. Model 2: Intelligent Re-tries 



Many variations are possible on the process that leads to Model 1. 

These are of interest because they represent possibilities for improved 
performance, or because they describe real system behavior. 

In the model of the present section we suppose that re-tries are in- 
telligent: a message in category R that attempts to seize a channel 

again does not do so at random, with a chance of destroying a message in 
progress. To represent this change, remove in (2.2) the possibility of the 
transitions involving a(t): (Q, R, L) ^ (Q+1 , R-1 , L), (Q-1 , R+1 , L), 

and ^ (Q, R-1, L+1 ) , replacing with the transition rate 



Q(t),R^-Q(t) + l,R(t)-l , probabil ity v(t)R(t)[l - -]dt (3.1 ) 

^ ^ ^ ~ c 



Effectively this change in Model 1 prevents losses, and, when a re-try 
prepares to access an occupied channel , it senses the presence of the on- 
going message, and refrains before destruction. 

This model can be analyzed by the technique described for Model 1, 
namely that of setting up approximating Ito-type stochastic differential 
equations, and then carrying out an expansion for c ^ ~ . The results are 
as follows. 

Deterministic Component : 

We find that for Model 2 

q'(t) = X(l-2q) - pq + vr(l-q) 
r' (t) = 2Aq - vr(l-q) 



(3.2) 



Stochastic Component : 

Use of the representation (2.4) together with the stochastic differ- 
ential equations and Ito's lemma yields 



i 

dX 


= 


f' \ 

X 


dY 


~t 


Y 


> 










(3.3) 
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the s.d.e. for a non-stationary Ornstein-Uhlenbeck process. In this case 
(suppressing t-dependence) , 






5t 



-(2X + y + 


vr) 


v(l-q) 




2X + vr 


-v^ 


-v(l-q) 

> 






-/yq 


/vr(l-q) 


0 


2Aq 


0 


-/vr(l -q) 



(3.4) 



(3.5) 



and 

Wt = (Wi(t), W^Ct))' . (3.6) 

a 4-dimensional standard Wiener process with independent components. 

Once again the equation (2.9) may be solved for the covariance function 

, with substituted for C^, and in place of . And once 

again we state the 

Result : (Q, R) is approximately bivariate normal as c ^ : 

(Q, R) ~ W(c(q,r), cE^) (3.7) 

Steady-State Results : 

If X(t) ->■ X , y(t) y , v(t) -»■ V as t ^ , then it may happen that 

q(t) -»■ q , and r(t) r , these satisfying the equations (3.2) with 

q' = r' =0 . The solutions are easily obtained, and are 



q 



r = = 2p\ 

v(y-x) n(l-p) 



(3.8) 



apparently it is necessary that q = p < 1 for steady state to occur. 

Note that in Model 1 it was necessary that q < 1/2 in order that steady 
state conditions occur. If we compare Models 1 and 2 when a = 1 we find 
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Model 1 



q = p 



(0 < p < 1) 



Mo del 2 

q = p , (0 < p < 1 ) 



r = 






221 



n(l-2p) " n(l-p) 

Clearly the setup of Model 2 results in considerably smaller expected wait 
per unit time (approximated by c(q+r)) than does that of Model 1. The 
improved service must be purchased in return for investment in the busy 
channel sensing capacity. The Model 2 can be compared to the model of Gaver 
and Lehoczky (1976). 

The steady state covariance matrix, Z , is obtained from (2.14), 
replacing C by A , and D by B . We find that 



0 

1 1 




b(a+b-d) 


2b^ 


b^ 




X + pq + vr(l-q) 


a 

1 2 




bd 


2ab 


ab 




-(2Xq + vr(l-q) ) 


a 

22 

L r' 




d^ 


2ad 


a(a+b)-bd 




4Xq + vr(l-q) 



(3.9) 



where 



a = 2A + p + vr , 
d = 2X + vr , 



b = v(l-q) 

A = 2ab(d-a-b) + 2b^d . 



(3.10) 
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4. Model 3: Transitory Version of Model 1. 



Consider next a transitory service system version of the basic Model 1. 
There are N messages to be transmitted, and each initially is transmitted 
independently and at a time having absolutely continuous distribution function 
F(t) , with F(0+) = 0 , and density f(t) . Once a given message is trans- 
mitted, no more appear from its particular source, so even though some 
messages enter the re-try population one or more times, the traffic gradu- 
ally fades away; the problem is fundamentally non-stationary . For analysis 
of a similar problem see Gaver, Lehoczky, and Perlas (1975). 

The following state variables are required. 



I(t) = the number of message arrivals that have occurred by time t ; 



0 < I(t) < N . 

Q(t) = the number of messages being transmitted at time t ; 

0 < Q(t) < c . 

R(t) = the number of messages in the re-try population at time t . 
We assume, as we did in formulating Model 1, that if a newly arriving message 
encounters a channel — selected at random — held by a message in progress, 
then both are instantly "destroyed," and immediately join R . No messages 
are ever lost. Messages in R re-try at rate v(t)R(t) . 

The transition scheme is given below. The individual message arrival 
rate at time t is seen to be X(t) = f(t)/[l -F(t)] . 



Transition 
(t to t+dt) 



Probability 



(I, Q, R) ^ (I+l, Q+1, R) 



Q 

X(t)[N - I(t)][l - ^]dt 



->■ (I+l, Q-1, R+2) 



Q 

X(t)[N - I(t)]^ dt 



^ (h 9 -''’ V 
^ (I. 9+1. R-I) 

^ (I, Q-1, R+1) 



y(t)Qdt 

Q 

v(t)R[l - 7]dt 

~9 

v(t)R^t . 



(4.1) 
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I 



From the above we can write down the approximating stochastic differential 
equations analogous to (2.3 ). These are, suppressing t-dependence state- 
ments , 

dl = X[N - I]dt + v/[N - I][l - |] dW, + /[N - Ip dW. (4.2) 

Q Q Q Q 

dQ = {A[N - I][l - - X[N - I]^ - uQ - vR^ + vR[l - -]}dt 



/ r / Q _ 

+ A[N - I][l - dW^ - /X[N - I]^ dW2 ■ AQ dW^ 



- /vR- dW, 
-c -4 



+ J\) 



R[1 - d«5 






Q Q Q 

dR = {2X[N - I]J + vR^ - vR[l- ^]}dt + 



2 



/ Q 

A^N - I] dW^ + 




Once again we study the behavior of the system as system parameters, in this 
case both N , the initial number of messages, and c , the number of 
channels become large. In fact, we relate these parameters as follows: 
c = 3N, 6 being a positive constant. Introduce the noise processes 
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I 



I(t) - Ni(t) 

X(t) = , 

M 

Q(t) - Nq(t) 

Y(t) = = . (4.3) 

R(t) - Nr(t) 

Z(t) = , 

where (i(t), q(t), r(t)) = lim (I(t)/c, Q(t)/c, R(t)/c) . 

N-x» 

Then an application of Ito's lemma and identification of terms of order 
and N yields the deterministic and stochastic components of the process. 

Deterministic Equations : 
i '(t) = X(1 - i) 

q'(t) = X(1 - i)(l - q/g) - X(1 - i)q/g - yq - vrq/g + vr(l - q/g) . 

(4.4) 

r'(t) = 2X(1 - i)q/g + vrq/g - vr(l - q/g) 

From the first equation and the definition of X it follows immediately that 
X(l-i) = f , the density of the arrival distribution. Note that when this 
substitution is made in the second and third equations the latter are precisely 
of the form of the corresponding equations of Model 1 : 

Correspondence Between Model 1 and Model 3 
Arrivals and Re-tries 

Model 1 Model 3 

X(t) f(t)/g 

v(t) v(t)/g 

In the present model f(“) = 0 , and hence the arrival rate is eventually 
zero, and some time afterwards the system is completely empty. The present 
model possesses no steady state solution, and to learn about system status at 
various times it is necessary to solve (4.4) numerically. 
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Stochastic Equations : 

The expansion technique shows that the noise process satisfies 






1 

dX 




X 


dY 


= 


Y 


dZ 




Z 



(4.4) 



where 



-X 

-X(l-2q/3) 

-2Xq/3 



and 



B. = 



/f(l-q/B) 

/f(l-q/B) 

0 



-vWF 



2/fq/B 



0 0 
-2f/B + y+ 2vr/B v(l-2q/B) 
2f/B + 2vr/B -v(l-2q/B) 

0 0 0 

-v4>rq/B /vr( 1 -q/B) 

0 



(4.5) 



/vrq/B -/vr(l-q/B) 



(4.6) 



and is a 5-dimensional standard Wiener process with independent 
components . 

Resul t : (I, Q, R) is approximately normal as N (hence c) ® : 

(I, Q, R) ~ W( c(f, q, r), cE^ ) 

Here satisfies the differential equation 






(4.7) 



and 



5t?; 



2fq/B 



f f(l-2q/B) 

f(l-2q/B) f(l-2q/B) + yq + vr -(2fq/B + vr) 

2fq/B -(2fq/B + vr) 4fq/B + vr 



(4.8) 
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As might be anticipated, the fact that 

Var[I(t)] = a (t) + F(t)[l-F(t)] (4.9) 

1 1 

may be deduced from (4.7). 
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5. Model 4: Transmission to Completion in Model 1. 

Suppose now, as may be quite realistic, that a message attempting to 
transmit on a channel does so until completion before it discovers that it 
has been "destroyed," i.e. garbled by others — which it also contributes to 
garbling. The analysis at once becomes much more complex because a channel 
may contain any number of destroyed messages. Our formulation is that of 
Model 1 , save for the change described above. 

The following state variables describe the system. 

Q(t) = the number of channels carrying good, i.e. ungarbled or unde- 
stroyed messages at time t . 

S|^(t) = the number of channels carrying exactly k messages that are 
destroyed at time t ; clearly S^(t) h Q(t) , and 

0^Q + ZS. < C. 

R(t) = the number of messages in R at time t . 

The transition probability scheme is summarized next. We write 

S+(t) = ? S. (t) . 
k=o'-'^ 
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Transitions 

(Q> ^1’ R) 

(Q+1 i S-j)..., R) 

(Q+1 . §1 , . . . , S|^, . . . , R-1 ) 

(Q~lj S-j)***) R) 

(Q“lj ^1’ ^2^^ >•••» R~l) 

(Q-1, Sp $2+1,..., S,^ R) 

(Q, S.J , . . . , $1^+1 , . . . , R-1 ) 
for k = 1 , 2, 3, ... 

^9’ 9l * ” " * ’ ’ ^k+l * * ■ ' * 

(Q, S-|~l,..., S|^,..., R+ 1 ) 

(Q , S-| , . . . , S|^+l , "1 , . . . , R+1 ) 



Probabilities 



Q+S. 

cX[l - ---- ] dt 
Q+S. 

vR[l - ^^]dt 
yQdt 

Q 

vR-^t 

Q 

S (5.1) 

rk 

vR — d t 
^k 

cX^t 

c 

-1 

y^t 

c 

^§k+i 



There will be a denumerable infinity of such transitions. While it is in 
principle possible to carry out the expansion technique, we shall content our- 
selves with a brief discussion of the deterministic equations analogous to 
(2.6). Following the example of Section 2, we can write down these differ- 
ential equations for the limits as c^ of Q/c , ^c , and Sj^/c (k=l,2,...), 
denoted by q , r , and : 
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q' = - (A + vr)0 - - q) - (p + A + vr)q 



r' = - vr + y J ks. 

k=1 (5.2) 

s,^ = y(k+l)s^^^ + (A+vr)S|^_.| - (k+A+vr)s,^ , k=l ,3,4,5, .. . 

$2 = Sps^ + (A+vr)(s^+q) - (2p+A+vr)s2 

The steady-state values, which exist under circumstances to be derived, are 
obtained by solving the above system of equations with the derivatives set 
equal to zero. To simplify, first divide through by y , putting 
p = A/y , n = v/y . After some algebra it is found that 

q = P (5.3) 

p + nr = 0 , so r = (6 - p)n'^ (5.4) 

where 6 is the smallest solution of the equation 

X "X tc c\ 

’Pi^ e = p , (5.5) 



provided one exists. 

S^ = p0 



e-V 

k! 



k > 2 



(5.6) 



Further analysis shows that in order for a steady-state to exist. 



p < e'^*^ ■ = 0.20588... 

+ 1 

Notice that this value is much smaller than the just-intolerable value of 
1/2 derived from Model 1. Clearly, transmission to completion of messages 
provides opportunity for many more transmissions to be destroyed. 
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6. Model 5: Transmission to Completion with Intelligent Re-tries 



A natural variation of the previous model is obtained by insisting 
on transmission to completion, but also allowing for an intelligent re-try 
capability, as in Model 2. This means that the transitions 

9l ’ 92''’^’"’ 5k’"’ ”^^9’ 5l’*"’ 5k+l’"*’ 

ruled out, i.e. have probability zero in this model. Consequently the tran- 
sition rates of (5.1) apply, with this change. Confining attention to the 
deterministic differential equations it may be shown that these have the 
form given below. 

q' = (X+vr)(l - - q) - (X+p)q 



r' = -vr(l - s - q) + y ks. 

k=l 

S|^ = (k+1) yS|^^^ + AS|^_^ - (A+ky)S|^ 

$2 = Sys^ + As-j - (A+2y)s2 + Aq . 

Once again let p=A/y, n=v/y, and solve for the steady-state values. 
These are 



( 6 . 1 ) 



q = p (6.2) 

, , I £m±£k!j:_ ( 6 . 3 ) 

l-(l+p)(eP-l) 

s^= p^ (6.4) 

k 

Sk= (1+p)^ ’ k=2,3,... (6.5) 

The condition for existence of a steady state is that the denominator of 
(6.3) be positive, which translates into the requirement that p < 0.50855... . 
The latter value may be contrasted to the just-tolerable value of unity 
derived for Model 2. 
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7. Numerical Results for Model 1. 



In order to check the quality of the diffusion approximations, a 
simulation program was written that realizes the two-dimensional Markov 
process describing Model 1. The latter was then exercised through 

c 

5 X 10 state changes for several values of the offered load, X/p , 
and for several channel numbers; the parameter a = 1 , and v/p = 0.08 
in this experiment. The results were generally supportive of the approxima- 
tion, as is seen by examining the following tables. Note, however, that 
assessments of the mean and variance of the number of active channels, Q , 
are generally in closer agreement than are those for the size, R , of the 
re-try population. 
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TABLE 1. 



Means and Variance-Covariance Values by 
Simulation and Diffusion Approximation 

C = 1 0 ; v/y = 2.0 





A/y: 0.48 


0.40 


0.32 


0.20 


E[Q] simulat. 


4.79 


3.99 


3.19 


1 .99 


di ffus . 


4.80 


4.0 


3.2 


2.0 


Var[Q] 


2.52 


2.52 


2.38 


1.76 




2.50 


2.52 


2.39 


1.77 


1 1 

1 1 

UJ 


65.13 


9.31 


3.29 


0.75 




57.6 


8.0 


2.84 


0.67 


Var[R] 


1470. 


50.26 


10.62 


1.57 




1466. 


42.81 


8.80 


1.35 


Cov[R,Q] 


0.034 


0.124 


0.166 


0.149 




0.038 


0.143 


0.177 


0.150 
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TABLE 2. 



E[Q] simulat. 
diffus . 

Var[Q] 

E[R] 

Var[R] 

Cov[R,Q] 



Means and Variance-Covariance Values by 
Simulation and Diffusion Approximation 



c = 25 ; v/y = 2.0 

X/\i: 0.48 0.40 0.36 



11.95 9.97 7.98 

12.0 10.0 8.0 

6.32 6.31 5.98 

6.26 6.30 5.97 

149.5 21.11 7.49 

144.0 20.0 7.11 



3524. 

3666. 

0.041 

0.040 



110.2 

107 

0.132 

0.143 



23.74 

22.0 

0.177 

0.177 



0.20 

4.99 

5.0 

4.41 

4.42 

1.73 

1.67 

3.56 

3.38 

0.150 

0.150 
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8. Summary and Conclusions 



In this paper it has been shown that an approximate approach, using 
stochastic differential equations, is effective for modeling an element of 
a complex communication system. The approximation improves as c , the 
number of channels available for transmission, becomes large. The adequacy 
of the model under such conditions is suggested by the theoretical results 
of Barbour (1974) and Kurtz (1971); numerical solutions of selected systems 
also attest to the adequacy of the approximation. The authors wish to 
gratefully acknowledge the research support of the Office of Naval Research. 
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